202201.00073v1 


chinaXiv 


摘 要 : 以 中 国 气象 局 提供 的 地 面 降 R 


ChinaXiv 合 作 期 刊 


tiwn 


ARID LAND GEOGRAPHY 


全 球 降雨 计划 GSMaP 5 IMERG 卫星 降雨 
产品 在 陕西 地 区 的 精度 评估 


FER, KR 昌 ”， KER 
(1. 陕西 省 地 表 系 统 与 环境 承载 力 重点 实验 室 ,陕西 西安 710127; 
2. 西北 大 学 城市 与 环境 学 院 地 表 系 统 与 灾害 研究 院 , 陕 西 西安 710127; 
3. 旱 区 生态 水 文 与 灾害 防治 国家 林业 和 草原 局 重点 实验 室 ,陕西 西安 710127) 


观测 资料 为 参考 ,利用 相关 系数 、 均 方 根 误差 和 相对 偏差 , 降 
雨 误 报 率 、 命 中 率 和 关键 成 功 指数 等 6 种 不 同 的 统计 分 析 指 标 ,从 年 、 季 、 月 、 日 4 种 不 同 的 时 间 尺 


度 对 两 种 高 分 辨 率 卫星 降雨 产品 (IMERG 和 GSMaP) 在 陕西 省 的 精度 进行 了 对 比 和 评价 ,并 对 二 者 
在 监测 强 降雨 过 程 中 的 表现 进行 了 对 比分 析 。 结 果 表 明 :(1) 在 年 尺度 上 ,GSMaP 的 数据 精度 高 于 
IMERG。GSMaP 和 站 点 观测 数据 高 度 相 关 ,而 IMERG 和 站 点 观测 数据 中 度 相 关 ;GSMaP 高 估 了 年 
尺度 的 降雨 ,而 JIMERG 低 估 了 年 尺度 的 降雨 (2) 在 季节 尺度 上 ,两 种 数据 均 在 夏季 的 精度 最 高 ， 


总 体 上 ,IMERG 在 季节 尺度 的 精度 亦 高 于 GSMaP。(3) 在 月 尺度 上 ,两 种 产品 与 地 面 观 测 数 据 均 呈 


现 较 高 的 相关 性 ,而 且 均 存在 一 定 程度 的 高 估 , 但 IMERG 比 GSMaP 更 具有 相对 较 高 的 精度 ,(4) 在 
日 尺度 上 ,GSMaP 的 数据 精度 略 高 于 IMERG。(5) 卫星 降雨 产品 的 数据 精度 与 降雨 量 有 关 ,总 体 表 


现 为 雨量 小 时 高 估 、 和 雨量 大 时 低估 。(6) LER 


产品 数据 精度 呈现 了 明显 的 地 域 差异 ,GSMaP 对 


陕西 省 的 降雨 总 体 表现 为 低估 ,其 中 对 陕 北 的 低估 最 为 明显 ;IMERG 对 关中 有 轻微 高 估 , 对 陕 北 和 
陕 南 则 存在 明显 低估 。(7) 通过 对 4 场次 强 降雨 事件 的 分 析 发 现 ,GSMaP 对 大 十 及 以 上 强 降 雨 事件 
的 监测 能 力 比 IMERG 略 强 。 研 究 结 果 可 以 为 该 地 区 的 气象 水 文 研究 在 选择 和 使 用 降雨 数据 资料 


时 提供 参考 。 
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文章 编号 : 


降雨 代表 了 地 球 表层 与 大 气 之 间 至 关 重 要 的 
能 量 和 水 分 交换 ' ,是 气象 气候 相关 人 研究 的 关键 变 
量 ”。 传 统 的 降雨 观测 方式 虽然 可 以 获得 小 范围 较 
为 准确 的 降雨 数据 ,但 是 存在 空间 分 布 不 均匀 以 及 
部 分 地 区 测 站 稀少 等 问题 ,难以 准确 反映 大 尺度 降 
雨 的 时 空 分 布 ””。 卫 星 遥 感 技术 能 够 克服 传统 观 
测 方法 在 时 空 上 的 局 限 性 ,实现 全 方位 .多 时 相 、 大 
范围 的 连续 降雨 观测 。1997 年 ,美国 发 射 了 热带 降 
十 测量 任务 5 (Tropical rainfall measuring mission, 
TRMM) E Æ ,其 搭载 的 降雨 雷达 (Precipitation ra- 
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dar,PR) 是 全 球 首 个 星 载 测 雨 雷达 ,可 以 进行 33"S~ 
35°N 之 间 的 降雨 观测 , 且 能 够 提供 暴雨 三 维 结构 ， 
这 对 精确 估计 降雨 起 到 了 重要 作用 “ ,为 大 尺度 水 
文 气象 的 研究 提供 了 新 途径 "”*"。 基 于 TRMM 提供 
的 多 年 降雨 观测 ,结合 其 他 人 研究 和 业务 卫星 提供 的 
降雨 产品 ,形成 覆盖 范围 为 50"S~502N 的 降雨 数据 ， 
即 TRMM 多 卫星 降雨 分 析 产 品 TMPA(TRMM multi- 
satellite precipitation analysis ) 9 ,其 空间 分 辩 率 为 
0.25° x 0.23"" 。2014 年 ,由 美国 航天 航空 局 (Na- 


tional Aeronautics and Space Administration, NASA ) 
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和 日 本 宇宙 航空 研究 开发 机 构 (Japan Aerospace Ex- 
ploration Agency,JAXA ) 联 合 实施 的 全 球 降 十 观测 
计划 (Global precipitation measurement, GPM) ,是 
TRMM 的 后 续 卫 星 降雨 计划 ”。 其 发 射 的 CPM 降雨 
卫星 携带 了 全 球 第 一 个 KuwKa 波段 双 频 测 两 雷达 
(Dual-frequency precipitation radar, DPR)""" ,能 够 更 加 
精准 地 监测 微量 降雨 (<0.05 mmh OMARE”, 
自 此 ,多 卫星 遥感 反 演 降雨 正式 由 TRMM 时 代 跨 入 
T GPM EHR GPM 时代 广 泛 应 用 的 两 种 卫星 降 
十 产品 为 IMERG (Integrated multi-satellite retrievals 
for GPM )™ 和 GSMaP (Global satellite mapping of pre- 
cipitation)" o + TRMM 48 HÉ , GPM 卫星 降雨 产品 履 
盖 范 围 更 广 (扩展 到 全 球 ) ,时 空 分 辨 率 更 高 (CS- 
MaP:1 h,0.1° x 0.1°; IMERG: 30 min,0.1° x 0.1°) "。 

近 些 年 来 ,国内 外 已 涌现 出 大 批 对 GPM 卫星 降 
雨 产 品 的 相关 研究。Beria 等 后 对 IMERG Final Run 
和 TRMM 3B43 在 印度 86 个 盆地 的 反 演 精度 进行 了 
对 比分 析 ,结果 显示 ,IMERG 在 整个 印度 盆地 的 所 
有 降雨 强度 上 的 表现 均 优 于 TRMM。Ning 等 “对 
GPM 产品 在 中 国 东 部 地 区 的 误差 特征 和 性 能 分 析 
中 表明 ,IMERG V04 对 中 国 大 部 分 地 区 的 降雨 有 显 
著 的 高 估 , 而 GSMap-gauged V06 对 降雨 的 估 测 偏差 
与 降雨 量 有 关 。 陈 汉 清 等 “在 对 GPM 卫星 降雨 产 
品 的 评估 工作 中 发 现 , 在 复杂 地 形 地 区 和 高 海拔 地 
区 ,GPM 卫星 降雨 产品 的 表现 不 理想 ,甚至 会 出 现 
不 可 徘 的 情况 , 且 其 精度 具有 明显 的 季节 差异 性 。 
王 思 梦 等 ”\ 李 媛 媛 等 "在 黑河 流域 和 黄河 流域 的 
研究 中 均 发 现 IMERG V04 存 在 对 高 海拔 地 区 降雨 
低估 的 问题 。GPM 已 应 用 于 流域 水 文 模拟 ,但 
是 ,卫星 降雨 数据 易 受 多 种 因素 影响 ,使 其 数据 精 
度 在 不 同时 期 不同 地 区 出 现 较 大 差异 。 因 此 ,在 
选择 降雨 资料 开展 区 域 性 研究 之 前 ,对 相关 卫星 降 
雨 产 品 的 精度 验证 十 分 重要 。 

陕西 省 地 跨 长 江 、 黄 河 两 大 水 系 ,地 形 气 候 复 
杂 。 万 相 均 等 2 根据 1981 一 2010 年 的 气象 站 点 实 
测 数据 ,分析 了 陕西 省 降雨 变化 的 时 空 分 布 ,发 现 
陕西 省 年 降雨 主要 呈现 南北 差异 型 的 空间 分 布 , 表 
现 为 陕西 南部 降雨 较 多 ALES BERN Be, (Ese 
对 陕西 秦 巴 山区 TRMM 3B42 的 精度 评价 结果 表明 ， 
TRMM 数 据 在 年 、 季 、 月 尺度 上 有 和 较 高 的 精度 ,日 尺 
度 精度 相对 较 差 ,降雨 事 件 探 测 能 力 表现 良好 。 
昭 昭 等 ”基于 地 理 加 权 回 归 (GWR ) 模 型 对 陕西 秦 


巴山 区 TRMM 3B43 数 据 进 行 降 尺 度 研 究 。 但 是 ,高 
ay HEE GPM 时代 的 卫星 降雨 产品 在 陕西 地 区 的 精 
度 还 不 清楚 。 理 解 从 TRMM 时 代 到 GPM 时 代 卫 星 
降雨 产品 在 空间 分 辩 率 及 覆盖 范围 提高 的 同时 ,其 
精度 有 何 变 化 "" ,对 于 进一步 促进 新 型 高 分 辩 率 卫 
星 降雨 产品 在 陕西 地 区 的 应 用 具有 重要 意义 。 

因此 ,本 研究 拟 对 两 种 CPM 时 代 的 卫星 降雨 产 
im (IMERG 和 CSMaP) 在 陕西 地 区 的 精度 进行 评估 
和 对 比 ,以 研究 区 内 地 面 观 测 数据 作为 参考 ,选取 
一 系列 精度 评价 指标 ,分 别 从 年 、. 季 月 尺度 和 日 尺 
度 对 两 种 卫星 降雨 产品 的 精度 进行 评价 ;同时 ,对 
两 种 数据 在 陕 北 、 关 中 和 陕 南 3 个 地 貌 类 型 区 的 精 
度 进行 横向 对 比分 析 ; 最 后 ,选取 典型 强 降 雨 事件 
分 析 两 种 数据 在 监测 强 降雨 事件 的 表现 。 


1 研究 区 概况 


陕西 省 位 于 中 国内 陆 腹地 ,黄河 中 游 ,邻接 8 省 ， 
具有 承 东 局 西 .连接 南北 的 区 位 之 便 , 介 于 105?29'~ 
111°15’E, 31°42'~39°35'N 之 间 , 由 南 向 北 依次 横路 
北 亚热带 、 暖 温带 和 中 温带 3 个 气候 带 ,南北 气候 差 
异 较 大 。 陕 西 降雨 总 体 呈 现 南 多 北 少 的 规律 六 , 受 
季风 和 地 形 影 响 ,陕西 的 降雨 季节 性 变化 特征 明 
显 , 总 体 表现 为 夏季 多 十 而 冬季 干燥 入 。 为 研究 卫 
星 降雨 产品 在 不 同 地 貌 类 型 区 的 误差 特征 ,以 北山 
和 秦岭 为 界 ,将 陕西 省 划分 成 陕 北 黄土 高 原 、 关 中 
平原 和 陕 南 秦 巴 山地 3 个 地 貌 区 (图 1)。 


2 数据 与 方法 


2.1 降雨 数据 

本 文 以 中 国 气 象 数 据 网 提供 的 月 值 数据 集 和 
日 值 数据 集 作 为 地 面 基准 数据 ,17 个 气象 站 点 的 空 
间 分 布 如 图 1 所 示 。 选 取 2014 年 12 月 一 2019 年 11 
月 的 数据 作为 年 尺度 分 析 的 数据 ,其 他 尺度 分 析 中 
使 用 的 数据 时 间 跨 度 均 为 2016 年 12 月 一 2017 年 
11 月 。 

本 文采 用 的 研究 数据 为 IMERG 和 GSMaP 两 种 
卫星 降雨 数据 。 根 据 处 理 过 程 的 不 同 ,IMERG 分 为 
Early、Late 和 Final 3 种 产品 。Final 产品 是 用 月 度 仪 
表 数 据 创建 的 研究 级 (Level 1~3) 数 据 , 其 精度 优 于 
Early 产品 和 Late 产 品 ”。 本 文 使 用 了 NASA 官 网 
提供 的 精度 较 高 的 IMERG Final, GSMaP 生产 了 
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图 例 得 到 

一 一 陕西 省 省 界 Ss eae 

si Bel PETIT 卫星 降雨 产品 的 格 点 数据 表示 该 格 网 范围 内 
o 信介 站 点 C N EX) PRAYER EPE I A Bt, E mm -d H Pa h ER 
高 程 /m x 5 ` 
mm 182-612 A mm +h? (小 时 产品 ); 地 面 站 点 记录 的 是 点 位 的 单位 
a 时 间 降 雨量 ,单位 也 是 mm.d-:! 或 mm.h ,所 以 两 者 
Gl 1271~2078 


E 2078~3594 


100 km 


注 :该 图 基于 国家 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 审 图 号 
为 GS(2019)3333 号 的 标准 地 图 制作 , 底 图 无 修改 。 下 同 。 
图 1 陕西 省 地 形 及 气象 站 点 分 布 图 
Fig. 1 Topography and distribution of meteorological stations 


in Shaanxi Province 


GSMaP_NRT, GSMaP_MVK 和 GSMaP_Guage 3 种 数 


在 单位 上 可 直接 匹配 。 在 空间 上 ,针对 地 面 站 点 的 
经 纬度 确定 其 所 对 应 的 格 网 ,由 于 所 用 站 点 较为 稀 
距 ,未 出 现 一 个 格 网 中 同时 存在 2 个 或 多 个 站 点 的 
情况 , 故 二 者 之 间 为 一 一 对 应 的 情况 。 文 中 所 使 用 
的 卫星 降雨 产品 数据 与 地 面 观 测 数据 均 由 数据 发 
布 方 进行 了 质量 控制 , 且 本 文 在 使 用 地 面 数据 之 前 
进一步 对 其 进行 了 检查 ,保证 了 数据 的 质量 。 
2.2 研究 方法 

文章 以 地 面 站 点 数据 为 基准 ,对 两 种 卫星 降雨 
产品 进行 站 点 检验 。 为 了 对 两 种 卫星 降雨 产品 进 
行 精 度 验 证 ,本 文选 取 相 关系 数 (Correlation coeffi- 
cient, CC) 、 均 方 根 误差 (Root mean squared error, 
RMSE) 和 相对 偏差 (Relative bias,BIAS)3 个 常用 指 
标 , 利 用 所 有 降雨 样本 对 卫星 降雨 产品 性 能 进行 定 
量 评价 。 在 年 尺度 ,基于 17 个 站 点 , 共 采 用 了 85 组 
样本 ;在 月 尺度 , 共 采 用 了 204 组 样本 ;在 日 尺度 , 则 
采用 了 479 组 样本 计算 这 3 个 精度 指标 。CC 可 以 衡 
量 卫 星 降 雨 产 品 估 测 值 与 地 面 观 测 值 之 间 的 相关 


据 产 品 ,CGSMaP_Guage 是 CSMaP_MVK sai 
GSMaP_MVK 的 基础 上 经 过 全 球 30000 多 个 CPC 
雨量 站 点 数据 校正 后 获得 的 1h 尺 度 的 cs 
mm 7), AAEH GSMaP_Guage 产品 ,数据 来 源 于 
JAXA 官 网 。 所 使 用 的 两 种 卫星 降雨 产品 的 参数 见 
#1. 


表 1 两 种 卫星 降雨 产品 主要 特征 参数 
Tab. 1 Important characteristic parameters of the two 


satellite precipitation products 
卫星 降雨 时 间 空间 (i T 数据 


产品 分 辩 率 由 ”分 辩 率 /(9) 范围 来 源 
IMERG 0.5 0.1 90°N~90°S_ NASA 
GSMaP 1.0 0.1 60°N~60°S JAXA 


Æ XXT GPM 的 两 种 卫星 降雨 产品 在 多 个 时 间 
尺度 上 进行 误差 评估 ,包括 :年 、 季 、 月 尺度 和 日 尺 
度 。 其 中 ,季节 尺度 的 数据 是 由 每 3 个 月 的 数据 累 
积 而 成 ,具体 为 :12 月 一 翌年 2 月 (冬季 ) 3—5 H 
(春季 ) .6 一 8 月 (夏季 ) .9 一 11 月 (秋季 )。 每 年 的 降 
十 数据 由 12 个 月 (12 月 一 翌年 11 月 ) 的 数据 累积 


程度 ; BIAS 用 于 衡量 卫星 降雨 产品 对 于 实测 降雨 值 
误差 的 平均 趋势 ;RMSE 反 映 了 卫星 降雨 产品 的 整 
体 误 差 水 平和 精度 。 其 中 ,CC 值 越 接 近 于 1, 其 他 2 
个 指标 的 值 越 接 近 于 0, 则 说 明 产品 误差 越 小 。 各 
项 指标 的 计算 公式 如 下 : 


YD,- DW,- W) 


CC=—— (1) 
[S057 So -my 
IIx 2 
RMSE = [Dy 2) 
2(W,-D)) 
BIAS = 三 一 x 100% (3) 
>>, 


式 中 ;为 样本 数量 ; D, 为 地 面 观 测 降雨 mm) ; D 
为 地 面 观测 降雨 平均 值 (mm); W 为 卫星 降雨 
(mm); W 为 卫星 降雨 平均 值 (mm )。 

分 类 统计 指数 可 以 衡量 卫星 降雨 产品 对 地 面 
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降雨 事件 发 生 的 监测 能 力 ,包括 降 雨 误 报 率 
(False alarm ratio, FAR) 命中 率 (Probability of de- 
tection, POD ) 和 关键 成 功 指数 (Critical success in- 
dex ,CSI)” ,其 中 ,POD 和 CSI 的 最 优 值 为 1,FAR 的 
最 优 值 为 0。 各 指数 的 计算 公式 如 下 : 


__F 

FAR HeT (4) 
_ H 

POD a (5) 
— H 

CST= (6) 


式 中 :代表 实际 未 发 生 而 卫星 观测 到 强 降雨 事件 
(卫星 估 测 值 z 阅 值 ,地 面 观测 值 < BB AR; H 
代表 强 降雨 事件 被 卫星 正确 观测 到 (地 面 观 测 值 > 


布 规律 。 与 由 站 点 观测 数据 插值 得 到 的 结果 (图 
2a) 相 比 ,IMERG 对 年 平均 降雨 量 明显 低估 ,而 GS- 
MaP 对 年 平均 降雨 量 明显 高 估 。 

在 年 尺度 上 ,2014 年 12 月 一 2019 年 11 月 
IMERG 和 GSMaP 两 种 产品 相对 于 气象 站 点 观测 数 
据 的 评估 指标 表明 ( 表 2) ,GSMaP 和 站 点 观测 数据 
高 度 相 关 (CC=0.89) ,而 IMERG 和 站 点 观测 数据 中 
度 相 关 (CC=0.56)。GSMaP 对 年 尺度 的 降雨 有 所 高 
估 (BIAS=13.60% ) , mi IMERG 对 年 尺度 的 降雨 有 所 
低估 (BIAS=-11.41%)。 综 上 ,两 种 产品 均 能 较 好 地 
反映 陕西 省 的 年 降雨 量 ,GSMaP 优 于 IMERG。 王 思 
梦 等 在 黑河 流域 的 研究 中 发 现 IMERG 在 年 尺度 
上 表现 为 低估 地 面 实测 降雨 ,这 与 本 文 结 论 相 背 。 


BIE, LÆ E) RZ MRR LEKA 
观测 到 而 实际 发 生 强 降雨 事件 (地 面 观测 值 > FL, 
卫星 佑 测 值 < 靖 值 ) 的 次 数 。 本 研究 以 25 mm: d'E 
AAA RA SE AP AY BB EL? ,利用 FAR、POD 和 CSI3 个 
指标 评价 和 对 比 两 种 卫星 降雨 产品 对 大 雨 及 以 上 
级 别 强 降雨 事件 的 监测 能 力 。 在 所 选 的 17 个 站 点 
中 共 选 取 了 66 组 满足 大 两 阀 值 条件 的 样本 计算 这 
些 精度 指标 。 


3 结果 与 分 析 


3.1 年 尺度 精度 评估 

由 基于 地 面 站 点 数据 Kriging 插值 和 基于 
IMERG CSMaP 得 到 的 陕西 省 2014 年 12 月 一 2019 
年 11 月 期 间 5 a 平均 的 年 降雨 量 空 间 分 布 ( 图 2) 可 
知 , 陕 西 省 的 年 平均 降雨 呈现 明显 的 南 多 北 少 的 分 


(a) 站 点 观测 数据 N (b) IMERG 数 据 

图 例 A 图 例 

Tee aimee eae 
-关中 与 陕 南 分 界线 ， 一 关 让 与 陕 南 分 界线 


p” 


- 260 


p” 


- 260 
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EE== 


这 可 能 是 由 于 黑河 流域 属于 干旱 区 ,降雨 量 较 小 。 
而 本 文 发 现 , 卫 星 降 雨 产品 的 数据 精度 与 降雨 量 有 
关 ,总体 表现 为 雨量 小 时 高 估 、 雨 量 大 时 低估 。 
此 ,对 于 降雨 量 较 小 的 黑河 流域 ,IMERG 对 其 降雨 
量 有 所 低估 。 
3.2 季节 尺度 精度 评估 

陕西 省 的 降雨 季节 性 变化 特征 明显 ,为 了 进 一 
步 分 析 IMERG 和 GSMaP 的 精度 是 和 否 在 季节 上 存在 
差异 ,有 必要 对 卫星 降雨 数据 在 季节 尺度 上 进行 精 
度 验证 ”。 本 文 分 别 统计 了 陕西 省 在 冬季 EE 
夏季 秋季 的 降雨 量 ( 图 3) ,并 比较 了 两 种 卫星 降雨 
产品 与 地 面 观测 值 的 CC、RMSE 和 BIAS 3 个 评估 指 
标 (图 4)。 

IMERG 和 GSMaP 的 拟 合 优 度 均 在 秋季 最 高 , 决 
定 系数 ( 尼 ) 均 达到 0.9 以 上 ,GSMaP 在 夏季 的 拟 合 优 


N (c) GSMaP 数 据 

人 图 例 

一 陕西 省 省 界 

一 陕 北 与 关中 分 界线 
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图 2 2014 一 2019 年 陕西 省 年 平均 降雨 量 空间 分 布 


Fig.2 Spatial distribution of annual precipitation in Shaanxi Province from 2014 to 2019 
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2 两 种 卫星 降雨 产品 相对 于 气象 站 点 观测 
数据 的 年 尺度 精度 评估 指标 
Tab.2 Annual scale accuracy evaluation indices of 
the two satellite precipitation products relative to 


meteorological station observation data 


卫星 降雨 相关 系数 均 方 根 误差 相对 偏差 

产品 (CC) (RMSE )/mm (BIAS)/% 
IMERG 0.56 77.69 -11.41 
GSMaP 0.89 92.65 13.60 


度 最 低 ( 尺 =0.3681) ,其 他 季节 的 拟 合 优 度 均 较 高 ， 
而 IMERG 在 冬季 的 拟 合 优 度 最 低 (R=0.2338) ,在 
夏季 的 拟 合 优 度 也 较 低 (R=0.4117)。 另 外 ,各 条 拟 
合 曲线 与 1:1 线 之 间 的 关系 (图 3) 表 明 , 卫 星 降雨 产 
品 在 不 同 季 节 相 对 于 观测 站 降雨 量 存 在 高 估 或 低 
佑 ,两 种 产品 在 春季 和 秋季 对 降雨 均 有 高 估 现 象 ， 
而 在 冬季 和 夏季 对 降雨 有 高 估 也 有 低估 ,这 表明 卫 


60 
(a) 冬季 y=0.8504x+12.476, 
£ 50 R?=0.7348 
1 40 
E 30 Z y=0.3066x+22.815 
ne R°=0.2338 
1 20 e IMERG 
w 10 e GSMaP 
Pi — 1:14% 
0 
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观测 站 降雨 量 /mm 
300 (c) BH y=0.4426x+213.04 
= 400 R=0.3681 。 
Pe e 
= : Rese 
xe ay Z y=0.5363x+163.2 
nF 200 R?=0.4117 
iin e IMERG 
到 100 。GSMaP 
0 一 1:1 线 
100 200 300 400 500 
观测 站 降雨 量 /mm 


星 降雨 产品 的 精度 可 能 与 降雨 量 的 大 小 关系 密切 ， 
降雨 强度 较 小 时 一 般 表现 为 高 估 ,在 降雨 强度 较 大 
时 则 多 为 低估 。 

IMERG 和 CSMaP 在 季节 尺度 上 均 表 现 出 较 高 
的 相关 性 (图 4) ,二 者 与 地 面 观测 数据 的 相关 性 在 
秋季 最 高 , 均 在 0.9 以 上 ;两 种 产品 的 BIAS 则 均 在 降 
雨 最 少 的 冬季 最 大 ,在 夏季 最 小 ,这 与 汪 梓 彤 等 " 
青藏 高 原 对 IMERG 的 多 尺度 精度 评价 中 得 到 结 
论 一 致 ,说 明 两 种 产品 均 对 夏季 降雨 的 高 低估 程度 
较 小 ,而 对 冬季 的 高 估 程 度 较 大 , 且 IMERG 在 夏季 
的 BIAS 为 负 值 ,说 明 IMERG 对 夏季 的 降雨 有 轻微 
的 低估 现象 ,而 CSMaP 的 BIAS 值 大 于 0, 说 明 GS- 
MaP 对 夏季 的 降雨 有 轻微 的 高 估 现 象 ; 两 种 产品 的 
RMSE 均 在 春季 最 大 ,在 冬季 最 小 ,说明 在 冬季 , 卫 
星 降雨 产品 数据 距 实 测 数据 最 稳定 。 两 种 卫星 降 


(b) Æ$ y=1.1628x+15.49 
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图 3 两 种 卫星 降雨 产品 季节 尺度 观测 降雨 与 地 面 实测 降雨 关系 


Fig. 3 Relationship between the observed precipitation of the two satellite precipitation products at 


seasonal scale and the measured precipitation at ground 
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注 :CC 为 相关 系数 ,RMSE 为 均 方 根 误差 ,BIAS 为 相对 偏差 。 下 同 。 
图 4 两 种 卫星 降雨 产品 在 陕西 省 内 季节 尺度 的 精度 评估 指标 


Fig. 4 Accuracy evaluation indices of the two satellite precipitation products at seasonal scale in Shaanxi Province 
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EENE: 全 球 降雨 计划 GSMap 与 IMERG 卫 星 降雨 产品 在 陕西 地 区 的 精度 评估 


十 产品 数据 精度 的 表现 不 一 ,可 能 与 陕西 省 降雨 的 
时 空 分 布 有 关 ,也 可 能 是 因为 不 同 季节 降雨 温度 和 
雷达 反射 率 有 不 同 的 变化 汪 ,使 卫星 降雨 产品 的 精 
度 表 现 出 明显 的 季节 性 差异 。 
3.3 月 尺度 精度 评估 

IMERG 和 GSMaP 两 种 卫星 降雨 产品 在 月 尺度 
上 相对 于 所 有 17 个 气象 站 点 观测 数据 的 评估 指标 


IMERG 在 降雨 量 小 于 100 mm 时 对 降雨 主要 呈现 高 
佑 状态 ,在 降雨 量 大 于 100 mm 时 ,对 降雨 主要 呈 低 
佑 状态 (图 5a) ;GSMaP 对 降雨 始终 呈现 高 估 状 态 
(图 5b)。 综 上 可 知 ,在 月 尺度 上 ,两 种 产品 与 站 点 
观测 数据 均 比 较 接近 ,IMERG 的 表现 相对 更 好 。 
3.4 日 尺度 精度 评估 

通过 雨季 (6 一 9 月 )IMERG 和 GSMaP 与 观测 站 


( 表 3) 表 明 ,两 种 产品 的 数据 和 站 点 观测 数据 均 呈 
现 出 较 高 的 相关 性 。CSMaP 的 均 方 根 误差 相对 较 
大 。 此 外 ,两 种 数据 的 BIAS 均 大 于 0, 反 映 了 在 月 
尺度 上 二 者 对 降雨 均 存 在 高 估 。 通 过 绘制 的 与 地 
面 观测 数据 的 散 点 图 (图 5$) 可 见 , 两 种 产品 的 估 测 
值 均 展 现 了 与 地 面 观测 值 较 强 的 一 致 性 。 此 外 ， 


表 3 两 种 卫星 降雨 产品 相对 于 气象 站 点 观测 
数据 的 月 尺度 精度 评估 指标 
Tab.3 Accuracy evaluation indices of the two satellite 
precipitation products relative to the observation data of 


meteorological stations at monthly scale 


卫星 降雨 产品 CC RMSE/mm BIAS/% 
IMERG 0.95 2.77 4.29 
GSMaP 0.94 7.41 11.46 
300 
g (a) IMERG 月 尺度 
# 250 
a 200 oor 
at Too. 
190 S y=0.9338x+7.0537 
z 100 R?=0.911 
= 50 --- 拟 合 曲线 
= 0 1:1 线 
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观测 站 降雨 量 /mm 


有 效 降 雨 事件 (日 降雨 量 >1 mm ) 的 散 点 图 (图 6), 结 
合 两 种 卫星 降雨 产品 的 评估 指标 ( 表 4) 可 以 看 出 ， 
GSMaP 和 地 面 观测 数据 的 相关 性 更 高 , 且 RMSE 和 
BIAS 都 更 小 ,由 此 可 知 ,日 尺度 降雨 中 GSMaP 的 数 
据 精 度 比 IMERG 的 高 。 两 种 产品 的 BIAS 均 呈现 负 
值 ,表明 在 日 尺度 上 二 者 对 地 面 降 雨 均 有 所 低估 。 
IMERG 和 GSMaP 均 对 大 于 40 mm 的 日 降雨 量 呈 现 
明显 的 低估 。 

按照 气象 规定 的 降雨 闵 值 ,将 降雨 事件 按 日 降 
雨量 分 为 小 十 (0.1 mm:d 必 降雨 量 <10 mm-d") 中 两 
(10 mm:d'< BERS <25 mm-d') AE (25mm-d'< 降 
雨量 <50 mm: d!) 大 十 以 上 (降雨 量 >50 mm-d"'), 
以 研究 不 同 降雨 强度 下 卫星 降雨 产品 的 误差 特 
征 。 两 种 产品 在 4 种 不 同 降雨 强度 下 的 评估 指数 


300 
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图 5 两 种 卫星 降雨 产品 月 尺度 观测 降雨 与 地 面 实 测 降雨 关系 


Fig. 5 Relationship between the observed precipitation of the two satellite precipitation products 


at monthly scale and the measured precipitation at ground 
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图 6 两 种 卫星 降雨 产品 日 尺度 观测 降雨 与 地 面 实 测 降雨 关系 


Fig.6 Relationship between the observed precipitation of the two satellite precipitation products 


at daily scale and the measured precipitation at ground 
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4 两 种 卫星 降雨 产品 降雨 量 观测 相对 于 气象 
站 点 观测 数据 的 日 尺度 精度 评估 指标 
Tab.4 Accuracy evaluation indices of the two satellite 
precipitation products to the measured data of 


meteorological stations at daily scale 


卫星 降雨 产品 CC RMSE/mm BIAS/% 
IMERG 0.50 3.75 -27.23 
GSMaP 0.62 2.38 -17.27 


(图 7) 表 明 ,两 种 产品 在 小 十 .中 十 .大雨 时 相关 性 
都 很 低 ,尤其 在 中 十 时 两 种 产品 的 相关 性 都 最 低 ， 
在 大 雨 以 上 降雨 时 相关 性 较 高 ;两 种 产品 的 RMSE 
均 在 小 雨 时 达到 最 小 值 ,同时 其 BIAS 均 大 于 0, 旦 
值 最 大 ,说明 两 种 产品 均 在 小 十 时 对 降雨 量 高 估 严 
重 , 并 且 GSMaP 对 降雨 量 的 高 估 程 度 比 IMERG 更 
大 一 些 , 而 二 者 的 BIAS 在 其 他 情况 下 均 小 于 0, 可 
以 看 出 两 种 产品 在 降雨 强度 较 大 时 对 降雨 存在 低 
fit , H. IMERG 的 低估 程度 更 大 。 以 上 结果 表明 ,在 
陕西 地 区 ,GSMaP 对 降雨 事件 和 降雨 量 的 捕捉 能 
比 IMERG 更 好 。 
3.5 不 同 地 区 精度 对 比分 析 

为 了 量化 日 尺度 下 两 种 卫星 降雨 产品 在 陕西 
省 不 同 地 貌 区 的 有 效 降 雨 的 误差 ,统计 了 IMERG 和 
GSMaP 两 种 产品 在 陕西 省 内 不 同 分 区 的 评估 指数 
(图 8) ,分 析 表 明 ,两 种 产品 在 关中 地 区 的 相关 性 均 
较 低 ,但 是 在 关中 地 区 的 RMSE 和 BIAS 均 最 小 。 


小 雨 ”中 雨 大雨 大雨 以 上 
雨 强 


IMERG 轻微 高 估 关 中 地 区 的 降雨 而 明显 低估 陕 北 
和 陕 南 地 区 的 降雨 ,而 CSMaP 不 同 程度 地 低估 了 这 
3 个 地 区 的 降雨 ,其 中 ,对 陕 北 地 区 的 降雨 低估 现象 
最 为 严重 (BIAS=-27.16% ) ,其 次 是 陕 南 地 区 (BI- 
AS=-13.99% ) ,对 关中 地 区 的 降雨 低估 程度 最 轻 
(BIAS=-1.48% ) 。 

以 上 分 析 表 明 ,在 日 尺度 降雨 中 ,卫星 降雨 产 
品 数据 精度 呈现 出 明显 的 地 域 差 异 ,GSMaP 对 陕西 
省 的 降雨 量 总 体 表现 为 低估 ,尤其 是 对 陕 北 地 区 的 
降雨 低估 最 为 明显 ;IMERG 则 是 对 陕 北 和 陕 南 地 区 
均 有 较为 明显 的 低估 。 曾 岁 康 等 、 王 思 梦 等 “在 
四 川 地 区 和 黑河 流域 的 研究 均 发 现 IMERG 对 高 海 
拔 地 区 降雨 表现 出 低估 现象 ,这 与 本 文 研究 结果 一 
致 。 而 两 种 卫星 降雨 产品 在 关中 地 区 的 表现 与 其 
他 地 区 的 表现 差异 较为 明显 ,一 方面 ,这 可 能 是 由 
地 域 导 致 的 气 修 地 貌 等 差异 导致 , 男 一 方面 ,不 同 
地 区 站 点 的 数量 可 能 也 会 造成 一 定 的 有 影响。 总 体 
来 看 ,GSMaP 的 数据 精度 高 于 IMERG。 
3.6 强 降雨 事件 观测 精度 对 比分 析 

为 了 对 比 两 种 卫星 降雨 产品 在 监测 强 降 雨 事 
件 时 的 表现 力 ,以 陕西 省 2017 年 4 场 暴雨 事件 ( 即 
2017 年 6 月 3 一 5 日 .7 月 25 一 27 日 .9 月 23 一 27 日 以 
及 9 月 30 日 一 10 月 4 日 ) 为 例 ,分 别 使 用 卫星 数据 和 
地 面 观 测 数据 进行 了 降雨 过 程 分 析 。4 场 暴雨 事件 
的 发 生地 主要 集中 在 陕 北 和 陕 南 地 区 ,2 个 地 区 在 
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图 7 两 种 卫星 降雨 产品 在 陕西 省 内 日 尺度 的 精度 评估 指标 


Fig.7 Accuracy evaluation indices of the two satellite precipitation products at daily scale in Shaanxi Province 
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图 8 两 种 卫星 降雨 产品 在 陕西 省 内 不 同 地 区 的 精度 评 佑 指标 


Fig. 8 Accuracy evaluation indices of two kinds of satellite precipitation products in different areas of Shaanxi Province 
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EUS: 全 球 降 十 计划 GSMaP 5 IMERG 卫星 降雨 产品 在 陕西 地 区 的 精度 评估 


所 统计 的 暴雨 事件 期 间 平 均 日 降雨 量 均 达到 60 mm 
以 上 。IMERG 统计 的 平均 日 降雨 量 为 24 mm, GS- 
Map 统 计 的 平均 日 降雨 量 为 37 mm, 两 种 产品 均 对 
降雨 量 有 所 低估 。 图 9 为 4 场 强 降雨 事件 逐 3 小 时 
降雨 变化 过 程 ,其 中 横 坐 标 均 表 示 第 n(n=0, 1, 2, …， 
24) 个 时 间 点 ,每 次 暴雨 事件 的 起 始 统计 时 间 为 凌 
晨 02:00 ,每 隔 3h 记 为 1 个 时 间 点 ,站 点 观测 的 降雨 
量 数据 是 由 日 降雨 量 除 以 24 得 到 的 平均 降雨 量 。 
IMERG 和 GSMaP 对 雄 雨 事件 过 程 的 描述 比较 一 致 ， 
暴雨 事件 出 现 峰值 节点 的 时 间 也 基本 相同 ,说 明 两 
种 产品 都 可 以 较 好 地 反映 暴雨 期 间 不 同 强 弱 的 降 
雨 阶 段 , 相 对 于 日 降雨 观测 数据 展现 了 更 多 的 降雨 
日 内 波动 细节 ,这 对 于 开展 更 细 时 间 尺 度 的 降雨 水 
文 过 程 研 究 具 有 重要 意义 。 不 过 两 种 产品 对 降雨 
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量 的 捕捉 仍 有 差异 ,尤其 是 在 峰值 处 ,二 者 的 估计 
量 相差 最 大 。 

为 了 比较 两 种 产品 对 强 降 雨 事件 的 总 体 捕捉 
能 力 , 以 25 mm-d' AEM, VEE A AR ee HE E DA 
上 的 所 有 降雨 事件 ,计算 POD FAR 和 CSI( 表 5)。 
结果 表明 ,IMERG 和 GSMaP 均 在 陕 北 地 区 表现 出 了 
低 命中 率 (POD ) 和 低 关键 成 功 指 数 (CSI) 的 情况 ,并 
且 CSMaP 在 陕 北 的 误 报 率 (FAR) 也 最 大 ,而 IMERC 
在 关中 的 误 报 率 最 大 。 男 外 ,IMERG 和 GSMaP 均 在 
关中 地 区 表现 出 高 命中 率 ,GSMaP 的 命中 率 在 关中 
地 区 达到 了 0.80, 误 报 率 也 最 低 , 只 有 0.20, 并 且 在 
关中 地 区 的 关键 成 功 指数 也 最 大 ,达到 了 0.67; 
IMERG 在 关中 地 区 的 误 报 率 最 小 ,为 0.33, 在 陕 南 
的 关键 成 功 指数 最 大 ,为 0.54。 据 此 推断 ,在 不 同 地 
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图 9 暴雨 期 间 降 雨量 变化 过 程 


Fig.9 Variation process of precipitation during rainstorm 


5 大 十 及 以 上 级 别 强 降雨 事件 中 两 种 卫星 降雨 产品 在 各 地 区 的 精度 评价 指标 
Tab. 5 Accuracy evaluation indices of the two satellite precipitation products in each region 


in heavy and above heavy rainfall events 


区 域 卫星 降雨 产品 命中 率 (POD) iets (FAR) 关键 成 功 指数 (CSI) 

BEAL IMERG 0.36 0.29 0.31 
GSMaP 0.55 0.54 0.33 

关中 IMERG 0.67 0.33 0.50 
GSMaP 0.80 0.20 0.67 

EETA IMERG 0.64 0.22 0.54 
GSMaP 0.72 0.35 0.52 
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Syd DX E Æ EN A m AGE E BE A H w 2 EE 
IMERG 和 GSMaP 均 在 关中 的 表现 更 好 。 

总 体 来 看 ,IMERG 和 GSMaP 在 陕 北 对 大 雨 及 以 
上 强 降雨 事件 的 监测 情况 最 差 ,对 关中 地 区 的 监测 
情况 最 好 ,并 且 GSMaP 对 大 十 及 以 上 强 降雨 事件 的 
监测 能 力 比 IMERG 强 。 


4 结论 


本 文 使 用 中 国 气 象 数 据 网 获取 的 时 间 跨 度 为 
2014 年 12 月 一 2019 年 11 月 的 月 值 数据 集 和 日 值 数 
据 集 作 为 地 面 基准 数据 ,利用 多 种 评价 指标 分 别 对 
IMERG 和 GSMaP 两 种 卫星 降雨 产品 在 年 . 季 月 尺 
度 和 日 尺度 上 进行 精度 评价 ,并 对 陕 北 、 关 中 和 陕 
南 3 个 区 域 进行 区 域 差 异性 分 析 ,得 到 如 下 结论 : 

(1) 在 年 尺度 上 ,GCSMaP 相 较 于 IMERG 表现 出 
更 高 的 精度 。GSMaP 总 体 表现 为 高 估 , 而 IMERG 总 
体 表 现 为 低估 。 

(2) 在 季节 尺度 上 ,IMERG 和 GSMaP 的 估计 精 
度 在 夏季 更 高 。GSMaP 在 不 同 季 节 与 地 面 观测 站 
数据 的 相关 性 都 比 IMERG 略 高 ,但 RMSE 和 BIAS 
偏 大 ,反映 其 估计 的 降雨 量 与 站 点 降雨 量 的 趋势 具 
有 一 致 性 ,但 波动 性 较 大 。 

(3) 在 月 尺度 上 ,IMERG 和 GSMaP 与 观测 数据 
的 相关 度 都 很 高 , 且 表 现 出 一 定 程度 的 高 估 ; 年 尺 
度 和 月 尺度 上 ,IMERG 均 比 GSMaP 对 降雨 量 的 估计 
更 好 。 

(4) 在 日 尺度 降雨 中 ,GSMaP 的 数据 精度 高 于 
IMERG ,但 其 数据 精度 均 呈 现 出 明显 的 地 域 差异 ， 
GSMaP 对 陕 北 地 区 的 降雨 呈 明 显 低 估 ,IMERG 对 陕 
北 和 陕 南 地 区 的 低估 程度 都 比较 严重 。 

(5) 对 不 同 降雨 强度 下 卫星 降雨 产品 的 误差 特 
征 的 研究 发 现 ,卫星 降雨 产品 的 数据 精度 与 降雨 量 
有 关 。 卫 星 降 十 产品 在 降雨 强度 较 小 时 一 般 表 现 
为 高 估 ,在 降雨 强度 较 大 时 则 多 为 低估 。 

(6) 在 陕西 地 区 ,GSMaP 对 降雨 事件 和 降雨 量 
的 捕捉 能 力 比 IMERG 好 ,IMERG 和 GSMaP 对 大 雨 
及 以 上 强 降 雨 事件 的 监测 能 力 在 陕 北 地 区 最 差 ,对 
关中 和 陕 南 地 区 的 监测 情况 相对 较 好 。 

(7) 通过 在 不 同 尺 度 对 两 种 卫星 降雨 产品 的 精 
度 进 行 对 比 和 评价 发 现 ,IMERG 和 GSMaP 均 对 陕 北 
黄土 高 原 区 和 陕 南 秦 巴 山地 区 的 降雨 量 存 在 明显 
的 低估 ,在 对 陕 北 地 区 强 降 雨 事件 的 估计 精度 


不 高 。 

卫星 降雨 产品 在 不 同 的 地 区 呈现 明显 的 精度 
异 ,其 精度 受到 气候 、 地 形 \ 传 感 带 、 降 雨 强度 等 
因素 的 影响 ”。 将 来 可 以 引入 地 形 .气候 等 因素 来 
对 产品 进行 算法 的 改进 和 误差 的 订正 "。 本 研究 
的 结果 可 以 为 陕西 省 的 气象 水 文 研究 在 选择 和 使 
用 降雨 数据 资料 时 提供 参考 。 
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Abstract: The Global Precipitation Measurement (GPM) mission is an international network of satellites that 
provide next-generation global observations of rain and snow. GPM was initiated by NASA and JAXA as a global 
successor to the Tropical Rainfall Measurement Mission (TRMM). In the GPM era, IMERG and GSMaP are two 
main satellite-based precipitation products. Compared with TRMM, GPM has wider coverage and higher spatiotemporal 
resolution. Understanding the accuracy and reliability of these products is important to further promote the 
application of new high-resolution satellite precipitation products. Therefore, this study intends to evaluate and compare 
the accuracy of IMERG and GSMaP in Shaanxi Province, China. In this paper, with surface precipitation 
observation datasets acquired from the China Meteorological Data Network used as reference, six evaluation 
indices, namely, correlation coefficient, root mean squared error, relative bias, false alarm ratio, probability of 
detection, and critical success index, were employed to evaluate and compare their accuracy at different temporal 
scales. Their performance in monitoring heavy rainfall events was also analyzed. Results show the following: (1) 
At an annual scale, the accuracy of GSMaP is higher than that of IMERG. GSMaP is highly correlated with the 
ground observation, while IMERG is moderately correlated with the ground observation. GSMaP overestimates 
the precipitation at an annual scale, whereas IMERG underestimates the precipitation. (2) At the seasonal scale, 
both data have the highest accuracy during summer. Overall, IMERG has a higher data accuracy than GSMaP. (3) 
At a monthly scale, the two products show a high correlation with the ground observation, but a certain degree of 
overestimation is observed. IMERG has a relatively higher accuracy than GSMaP. (4) At a daily scale, the 
accuracy of GSMaP is higher than that of IMERG. (5) The accuracy of satellite precipitation products is related to 
total precipitation volume, which is generally overestimated when the precipitation volume is small and 
underestimated when the volume is large. (6) The accuracy of satellite precipitation product data shows obvious 
regional differences. GSMaP underestimates the precipitation in Shaanxi Province as a whole, especially in 
northern Shaanxi. IMERG slightly overestimates the precipitation in Guanzhong, but significantly underestimates 
the precipitation in northern and southern Shaanxi. (7) An analysis of four heavy precipitation events shows that 
the monitoring capability for heavy precipitation events of GSMap is slightly stronger than that of IMERG. This 
study is expected to provide a reference for the selection and application of precipitation data in meteorological 
and hydrological studies in this area. 
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